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1. Introduction The approaching date for broad band gravity wave detectors to come 
online makes the study of astrophysical sources of gravitational radiation immediately press- 
ing. Collisions among black holes and neutron stars are especially likely to generate detectable 
gravity wave signals. The binary black hole collision problem which suggests itself as the nat- 
ural "two-body problem of general relativity" is particularly attractive to a theorist since it 
allows us to restrict consideration to vacuum spacetimes. Much progress has been made to- 
ward solving this problem numerically in recent years, but the goal of a complete calculation, 
encompassing the inspiral of the holes, merger and the final ringdown of the remaining black 
7— I ■ hole still eludes us. Numerical techniques have been successfully applied to calculating the 

final stages of the collision, but there are strong limitations on the evolution time of the nu- 
merical codes as well as on the numerical domain size that these calculations can support. > 
This means that these calculations must begin from a configuration in which the black holes 
^C) ' are already close. 

To apply general relativity to astrophysical phenomenology it is natural invoke the 3+1 
formulation. Then phenomenological problems divide neatly into two steps. We must first 
\Q • define initial data on a spacelike Cauchy slice which satisfies Einstein's equations on that slice 

and which is appropriate for the physical problem under consideration. Then we must use 
the remainder on Einstein's equations to integrate that data through time. Once the initial 
data have been specified the procedure to calculate the evolution is well defined in principle 
even if it may be quite difficult in practice. Thus, in principle, selection of the initial data 
] specifies the entire spacetime. For problems of astrophysical interest it is then paramount to 

begin with initial data which correspond to a realistic astrophysical configuration. If the black 
q-( holes begin very far apart we could approximate each as a stationary spherically symmetric 

black hole and there would be no ambiguity in the initial data. But once the black holes are 
close enough to interact significantly it is difficult to say what are astrophysically realistic 
initial data. Most of the numerical calculations for binary black hole collisions that have 
. £h ! been attempted thus far have begun evolving from initial data constructed via the conformal 

flatness formalism. A special case of this, the so called Bowen-York method 3 specifies how 
to derive a family of solutions to the initial value equations which have a conformally flat 
spatial metric and vanishing mean curvature and which can be interpreted as boosted or 
spinning black holes. There are problems with the assumption of conformal flatness though, 
particularly for spinning black holes. For an initial data set to be astrophysically realistic, 
a minimum requirement is that sufficiently separated black holes look like Kerr black holes. 
Since we do not know of a foliation of the Kerr spacetime which admits conformally flat 
spacelike slices, we are unable to exactly represent Kerr black holes in such a formalism. This 
difference between Kerr black holes and Bowen-York spinning black holes is demonstrated by 
the nonzero gravitational radiation emitted by a single black hole. 4 

In this paper we describe a novel method for solving the initial value equations of general 
relativity for a broad class of 3-metrics which includes axisymmetric slices. Unlike the Bowen- 
York approach, where the extrinsic curvature is selected first and then the metric is calculated 
from the constraints, we first choose the metric and then calculate the extrinsic curvature. 
After introducing our technique, we apply it to the Kerr black hole demonstrating that the 
Kerr hole does indeed come out as a special case. By comparison with this solution, we 
prescribe a method which may be used to introduce multiple Kerr-like black holes on the 
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initial hypersurface. In particular, if the black holes are sufficiently separated, then in the 
neighborhood of each one, the data approach those for a Kerr hole. Moreover, if the two holes 
are sufficiently close, the data also resemble that for a single black hole. 

2. Solution of the Constraints In this section, we shall find the general solution to the 
constraint equations for a Cauchy slice whose intrinsic geometry admits a surface-orthogonal 
Killing field. In particular, this is true of the surfaces of constant Boyer-Lindquist time in 
the Kerr spacetime. Later, we will construct initial data with spatial geometry satisfying this 
same condition for two black holes. 

For definiteness we will assume that our Killing vector field corresponds to an axial 
symmetry, though this restriction on the topology of the is not essential to our construction. 
Let <j> be the azimuthal coordinate so that d/dcf) is the Killing field. The existence of orthogonal 
surfaces implies that our spatial metric can be cast in the form q i jdx t dx : > — hijdx t dx i + B d(p 2 , 
where the two-metric hij is restricted to surfaces of constant <f>. Since these surfaces have the 
topology of half-planes, we can write hij = Ahij where hij is a flat metric and A is a suitable 
conformal factor. If we choose two of the three spatial coordinates to be constant along 
the orbits of the isometry, then the components of hij in our coordinate system as well as 
the functions A and B will not depend on <p. If we choose rectangular coordinates in the 
half-plane, then our metric assumes the form 

qij dx i dx j = A(x, y){dx 2 + dy 2 ) + B(x, y)d(f> 2 . (1) 

The initial value equations of geometrodynamics are a set of four coupled partial dif- 
ferential equations which relate the metric and extrinsic curvature Kij of a hypersurface 
embedded in a vacuum spacetime. In tensor notation they are written as a vector equa- 
tion known as the momentum constraint and a scalar equation known as the Hamiltonian 
constraint: 

V j K ij =d i (q jk K jk ) 

{q ac qbd _ q ab q cd )KabKcd = (3) n[q] 

From here we will proceed as follows. First we will show that the momentum constraint 
reduces to a set of Cauchy- Riemann equations to be solved for K xx , K yy , and K xy and the 
condition that K X( p and K y ^, can be derived from a potential. Then the Hamiltonian constraint 
provides a condition which determines this potential. In the coordinate system introduced 
last paragraph, the momentum constraint looks as follows: 

dJBiK^ d(BiK xy ) = OjBhq _ IdQogjB/A)) H 

dx dy dx 2 dx { xx vv> 

d(BhK xy ) d{Bh Kyv ) djBiK) l dQog(B/A)) i 

dx + dy =A ^y—-2 d~y Ba ( K ** + K w) 

d(BhK x4> ) | djBiKyj) _ Q 
dx dy 

The general solution for the third equation can be written at once in terms of a potential 
function u(x, y): 

K x4> = +B-i^ 

Q y (2) 

The other two equations involve the remaining four components of the extrinsic curvature. 
Since there are two more quantities than equations, we may take two of the quantities as 
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known and solve for the remaining two. It proves advantageous to take the quantities P and 
Q defined as P := B?(K XX — K yy )/2 and Q := B^K xy as our unknowns since our equations 
then take the following form: 

( dp dQ 

dx- + dy- =j{x > y) 

dQ dp 

where the sources j and k are defined as 

d(BiK) l d(log(B/A)) i \ d{B\K xx + BhK vv ) 

J ' 5a; 2 5a; X:E wJ 2 dx 

d(B^K) l dQog{B/A)) i + (3) 

fc = " 2 5 (Jfra + + 2 

A" := i4(JT X!t + Kyy) + BK H 

These equations have the form of Cauchy-Riemann equations with source and can hence 
be solved by well-known methods of potential theory. For instance, we might write P and Q 
as integrals of j and k against Green's functions. 

P(x,y)=f ds {G n {x,y,x',y')n(x',y') - Gi(x,y,x',y')m(x',y')} 

JdD 

dx'dy' {G n {x,y,x',y')k(x',y') - Gi(x,y,x' ,y')j(x' ,y')} 

JD 



Q{x,y)=f ds {GK,(x,y,x\y')m(x\y') + Gi(x,y,x' ,y')n(x' ,y')} 

JdD 

dx'dy' {G R (x, y, x' , y')j(x', y 1 ) - Gi(x, y, x , y')k(x , y')} 

/d 

K xy = B~iQ 



L 



(4) 



K XX = ^J^ + B^P 

The choice of Green function will depend on the domain D over which x and y range and 
the conditions imposed on its boundary. In addition to this solution of the inhomogeneous 
equations, we might also require a solution of the homogeneous system. For applications 
related to the two-body problem, it is natural to consider the half-plane y > as our domain 
and ask that the extrinsic curvature go to zero as y — > oo. With these specifications, the 
relevant choice of Green function is given as 



G R (x,y,x',y') 
Gi(x,y,x',y') 



y'[(x - x') 2 + y 2 - x 2 



n[(x - x' f + (y - y') 2 ][(x - x') 2 + (y + y') 2 } 

lyy'jx - x') 

n[(x - x') 2 + (y- y') 2 ][{x - x') 2 + (y + y') 2 ] 



Before proceeding further with our calculations, let us pause to consider our ultimate goal. 
We began by specifying a metric by a choice of A and B. Next, we chose K and K xx + K yy . 
The solution of the equations above for P and Q determine two more components of the 



3 



extrinsic curvature. The only thing we still need for a complete specification of the initial 
data is the u potential. We shall now proceed to determine u by solving the remaining initial 
value equation, which is the Hamiltonian constraint. 

The Hamiltonian constraint can be regarded as an equation for the hitherto arbitrary 
function u(x,y). Substituting formula 2 into the Hamiltonian constraint and rearranging, we 
arrive at the formula 

' du\ /<9u N 
dx) \d V/ 



= 2v{x,y) (5) 



9 d2 

2v(x, y) := AB 2 - — [{K xx f + K xx K yy + (K yy ) 2 + (K xy f 

(3) 1 J>£?* d 2 Bi\ 1 f d 2 (\ogA) d 2 (logA) 

2ABh \ dx 2 dy 2 J AA\ dx 2 dy 2 

Equation 5 has the form of a Hamilton-Jacobi equation for a particle moving in a two- 
dimensional potential. As is well known from analytical dynamics, such a partial differential 
equation can be solved using the method of characteristics. Given a solution u(x, y) to the 
equation, the characteristics are a congruence of curves with the property that the tangent 
vector at any point on one of these curves equals the gradient of u(x,y). If we label the 
various curves in the congruence by a parameter s and describe them by the equations 

j x = X(t;s) 
\y = Y(t;s) 

then this property is expressed analytically as 

du _ dX 

dx dt 
du_dY_ 

dy dt 

From these equations above and the non-linear partial differential equation 5, it follows that 

d 2 X dv 



dt 2 dx 
d 2 Y dv 
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dt 2 dy 

To solve the partial differential equation, we can reverse the logic of the last paragraph. 
Start with a one parameter family of solutions to the second order ordinary differential equa- 
tions 6. (These solutions can be thought of as the trajectories of a fluid of particles moving in a 
potential well —v(x, y).) Differentiating with respect to t, obtain a vector field (dX/dt, dY/dt). 
This vector field will have vanishing curl and, by integrating, one can obtain the scalar u(x, y) 
of which it is the gradient. For our purposes, however, only the gradient of u and not u itself 
is required, so there is no reason to actually carry out the integration. Recalling equation 2, 
we can write the following formula for two components of the extrinsic curvature tensor: 

dt 

K B-h dX 
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Along with formulas 3 and 4, the above formula can be used to compute the extrinsic curvature 
given the four freely specifiable functions A(x, y), B(x, y), K(x, y), and K^(x, y) and having 
obtained a one-parameter family of solutions to the system 6. 

3. Single black Hole In the last section, we described a method for solving for the ex- 
trinsic curvature that involved a Hamilton-Jacobi equation and mentioned that the method is 
applicable to Kerr spacetimes. In this section, we will consider the geometry of a hypersurface 
in this spacetime from the standpoint of our method. In the next section, we will utilize this 
representation of the initial value problem for a Kerr spacetime to construct initial data for 
the two black hole problem. 

A stationary and axisymmetric spacetime admits a preferred foliation in which time is 
given by locally non-rotating vector field. For the case of Kerr, this is known as Boyer-Linquist 
time. In terms of our formalism the intrinsic geometry on such a constant time slice of Kerr 
can be written as 

qij dx l dx j = Ahij dx l dx 3 + B d4> 2 

r 2 + a 2 cos 2 9 



A = 



R 2 

r 2 + a 2 2mra 2 sin 4 9 



B = — + Rj ™ 

hij dx'dx 3 = dR 2 + R 2 d9 2 

R := ^ r - 2m + 2\/m 2 ~ - a 2 + \J r - 2m - 2\J m 2 

where r is the usual Boyer-Lindquist radial function and R, 9 are polar coordinates in the 
half-planes of constant <f>. The metric hij in these planes is seen to be manifestly flat. 

Since we will be considering multiple black holes, it may be helpful to discuss this form 
for the spatial part of the Kerr metric a little further. If A = B, the spatial metric would be 
conformally flat. That A and B differ makes the conformal flatness formalism inapplicable. 
It is interesting to note the limiting behavior of these functions. 

A=l+ 2 f + 0{^) r^oo 

B > A everywhere 

§ -» 1 + O (a 2 ) a -» 

| -» 1 + O (sin 2 6>) 9 -> 

Thus we see that B =/= A, but differs from A minimally except in a region near each hole 
and off the symmetry axis. Note that this description of the Kerr data is similar to that 
studied by Brandt and Seidel 6 . They studied initial data where the Kerr spatial metric was 
modified by the addition of a Brill wave, which changes the ratio Bj A. 

In Boyer-Lindquist coordinates, the only non- vanishing components of the extrinsic cur- 
vature for a hypersurface of constant t are K r $ and Kq^. From formula 2 of last section, 
we know that these can be expressed in terms of a potential function u(r, 6). Doing the 
computation, we find that the potential responsible for these two non-vanishing components 
is 

u(r, 9) = ( / + C 9 ° S %- (r 2 cos 9 - a 2 cos 9 - 2r 2 ) 
r 2 +a 2 cos 2 9 v ' 

We know that this potential function is the solution to a Hamilton-Jacobi equation. To 

find out which solution, we can inquire into the congruence of characteristic curves with which 
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it is associated. Recall that these curves are orthogonal to the level curves of u. Having the 
explicit expressions at hand, we could formulate this orthogonal trajectory problem concretely 
and attempt to solve it. However, for our purposes such details are unnecessary, and we shall 
content ourselves with a qualitative discussion. 

First, let us note that in the region r > a and away from the axis of symmetry, u is 
non-singular and its gradient does not vanish. As a consequence, the curves we are interested 
in will be smooth and free from bifurcations and fixed points. To determine their topology, we 
shall consider the rays 9 = 0, ir/2, n. On the ray 9 = n/2, u is constant, so the characteristic 
curves must be orthogonal to this ray. Examining the metric h, we see this means that they 
must point along the 9 direction. Near the symmetry axis, 9 = 0, the gradient of u goes to 
zero, so we shall make an expansion in powers of 9. 

V 2 - n 2 

u(r,9) = -4 +I - TTI -,9* + 0(9Z) 

From this, we see that the gradient of u is parallel to the symmetry axis for small 9 and that 
its norm vanishes like 9 3 . Hence the orthogonal trajectories will point at right angles towards 
or away from this axis. By symmetry, the same situation will occur near the other half of 
the axis, 9 = tt. Combining these observations, we conclude that the various curves of our 
congruence start out perpendicular to one end of the symmetry axis, arch around the black 
hole, and return perpendicular to the other end of the symmetry axis. 

To generate this initial data set for Kerr metric within our formalism, we first have 
to specify the spatial metric by choosing A and B as given in equation 7. Then, making 
the choices K = = and specifying that the congruence of curves needed for the 

construction of K r( p and Kg^ meet the symmetry axis perpendicularly will guarantee that the 
correct extrinsic curvature comes out of the construction. The one remaining choice is the 
overall sign of the extrinsic curvature, which determines the sign of the angular momentum. 
This choice can be made by specifying on which side of the symmetry axis the characteristics 
originate. 

4. Two Black Holes Given an axially symmetric spatial metric, our formalism shows 
how to find an extrinsic curvature which will satisfy the initial value problem of general 
relativity. Our strategy for representing two rotating black holes will be to choose a spatial 
metric which resembles the spatial metric of the Kerr solution in three regions. These are the 
two regions near the chosen locations for the black holes on the symmetry axis and the region 
near spatial infinity. Then we will use the methods described above to calculate an extrinsic 
curvature which satisfies the initial value problem. In general, there will be many solutions 
from which we will select one that is consistent with the description of the Kerr extrinsic 
curvature whenever the spatial metric approaches that of the Kerr solution. 

The first step is to construct a three-metric which has the appropriate behavior in the 
near and far limits. One way to do this is to use the form given above for A and B of the 
Kerr metric and replace ^, ^, and sin0, with functions which approach the correct values 
in the limits near one of the holes or away from both. Let the Kerr metrics approached near 
each of the two holes have parameters {mi, a\) and (777,2, 02). Then we want the parameters 
of the metric approached at infinity to correspond to an appropriate sum of the masses and 
angular moments associated with the two holes. Thus we want 

TO far = m i+ m 2 
and 



a far m far = a i m i + a 2^2 
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More explicitly then, if we rewrite the quantities A and B for the Kerr metric with m.f ar 
scaled to unity and replace R, a, and sin# with R, a, and s, which are functions of x and y, 
we have 

r{R) 2 + a 2 {l~~s 2 ) 



A = 



B 



R 2 

'r{R) 2 +a 2 
R 2 



+ 



2rR + a 2 s 2 



R A 



2 • 2 

r sin i 



We leave the expression U R 2 sin 2 0" appearing in B unchanged since this gives the distance 
from the symmetry axis and does not depend on the location of the black holes. In order for 
our initial data to behave properly in the three limits, we specify the following limits on the 
new functions: 



mi 



R 



s 2 



R -» 


JL, 




mi 


S 2 -> 


fl 2 sin 2 9 






a — ► 




mi 


R -» 


ik 


m 2 


s 2 -> 


R 2 sin 2 S 




J?2 



Q2 
m 2 



sin 



aimi+aim2 




m 2 



and 



oo 



mi 



m 2 



Here i? and 9 are polar coordinates on the conformally flat two-space and R\ and R2 are 
the distances from the chosen locations of the black holes in the flat metric. The third limit is 
the same as the first with 1 and 2 interchanged. A natural choice for the function R(Ri, R2) 
with suitable limiting behavior is 



R = 



mi 1 in 2 
Ri " r R 2 



It can be seen from equation 7 that in the case of non-spinning black holes, where we expect a 
to vanish identically, B becomes equal to A and A reduces to the conformal factor associated 
with the familiar Brill-Lindquist initial data for non-spinning black holes. 
We can satisfy the conditions for s 2 by choosing 

_ 2 R 2 sin 2 

s 2 = = 

R 2 

In addition to the requirements given above, there is a much less obvious condition that 
we must impose to assure that we can specify a Kij which is consistent with the Kerr extrinsic 
curvature. In the case where the angular momenta for the two black holes point in opposite 
directions, we will need to ensure the existence of a surface separating the two black holes on 
which the extrinsic curvature vanishes. This is essentially because needs to change sign 
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in passing from the region near one black hole to the region near the other. If the extrinsic 
curvature vanishes, the Hamiltonian constraint forces the scalar curvature to vanish as well. 

A straightforward way to guarantee the existence of such a surface on which the curvature 
scalar vanishes is to require that the gradient of a vanishes whenever a vanishes. This follows 
from the following observations: 

1) With our choice of R, the scalar curvature vanishes when a vanishes identically. 

2) The metric depends on a only in the form a 2 . The only terms surviving in the expression 
for the scalar curvature evaluated at a point where a = have the factor q^aihj. 

A function a which meets our requirements is 



These choices for R, s, and a specify a metric which meets our requirements and we will 
take this as our spatial metric for the rest of the discussion. There are certainly other equally 
well-qualified and well-motivated choices of spatial geometry satisfying the conditions we have 
given. 

Having specified a three-metric, we now want to find an extrinsic curvature which solves 
the constraints and resembles the extrinsic curvature for the Kerr case in the appropriate 
limits. Since K zz and the mean curvature K vanish in the Kerr case, we again choose them to 
vanish. This means we use the trivial solution to the equations 3 and 4. Then what remains 
is to find a solution of equation 5 which reduces to the Kij for Kerr in the limits. 

To determine K rt p and Kg^, specify that, as in the Kerr case, the characteristics meet the 
symmetry axis perpendicularly. Then the choice of sign of the tangent vector field along these 
curves must be made consistent with the desired sign for the angular momenta at each hole. 
In the case that a\ and a2 differ in sign, there are two regions separated by a curve on which 
a vanishes. Since we have demanded that ^TZ also vanish here, we can choose the signs of 
the vector field independently in each of these regions without introducing a discontinuity in 
the extrinsic curvature. 

5. Conclusions 

Here we have described a method for generating initial data on a Cauchy slice with any 
3-metric that admits a surface-orthogonal Killing field. Given such a metric, our method 
generates all extrinsic curvature tensors which are invariant under the motion generated by 
the Killing field. In particular, this gives an affirmative solution to the "half-sandwich" 
conjecture for this class of spatial metrics and perhaps our method will find use in studying 
the "thin-sandwich" problem. 

We have also described a set of solutions to the initial value problem which can be 
interpreted as two black holes arranged along a common axis of symmetry. Although we 
have discussed only the two black hole solutions, it should not be difficult to generalize the 
technique to an arbitrary number of black holes. 

Our solution has the advantage over the Bowen-York solution that the geometry in the 
vicinity of each singularity does approach the Kerr geometry when the black holes are well 
separated. We suggest these data as alternatives to the Bowen-York data for exploration 
because they represent well-separated spinning black holes. Still, for the case of nearby black 
holes, our choice of initial conditions includes a certain amount of arbitrariness that can only 
be resolved by a more careful consideration which takes into account the precise evolution 
which has brought the black holes close to each other. 
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Abstract 



This work consists of two distinct parts. In the first part we present a 
new method for solving the initial value problem of general relativity. Given 
any spatial metric with a surface orthogonal Killing field and two freely spec- 
ified components of the extrinsic curvature we solve for extrinsic curvature's 
remaining components. For the second part, after noting that initial data for 
the Kerr spacetime can be derived within our formalism we construct data 
for axisymmetric configurations of spinning black holes. Though our method 
is limited to axisymmetry, it offers an advantage over the Bowen-York pro- 
ceedure that our data approach those for Kerr holes in the limit of large 
separations and in the close limit. 
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